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We analyze a model of resonant point-contact tunneling between multiple Luttinger 
liquid leads. The model is a variant of the multi-channel Kondo model and can be 
related to the quantum Brownian motion of a particle on lattices with 7r-fiux through 
each plaquette (in the 3-lead case, it is a honeycomb lattice with 7r-flux). By compar- 
ing the perturbative and instanton gas expansions, we find a duality property of the 
model. At the boundary, this duality exchanges Neumann and Dirichlet boundary 
conditions on the Tomonaga-Luttinger bosons which describe the leads; in the bulk, 
it exchanges the 'momentum' and 'winding' modes of these bosons. Over a certain 
range of Luttinger liquid parameter, g, a novel non-trivial intermediate coupling fixed 
point controls the low-energy physics. The finite conductance at this fixed point can 
be exactly computed for two special values of g. For larger values of g, there is a 
stable fixed point at strong coupling which has enhanced conductance resulting from 
an analogue of Andreev reflection at the point contact. 

PACS: 73.23.-b, 71.10.Pm, 73.40.Hm, 73.40.Gk 



I. INTRODUCTION 

Despite being a subject of intense interest in recent years, the study of strongly correlated electron systems has had 
a checkered history, primarily for two reasons. On the one hand, non-perturbative techniques - of which there are 
precious few - are required for their analysis. At the same time, these systems often exhibit unexpected phenomena, 
rendering useless our intuition culled from Fermi liquid theory and other essentially perturbative problems. The 
greatest progress has been made on one-dimensional systems and, particularly, quantum impurity problems. In this 
arena, powerful techniques such as conformal field theory Q and the Bethe ansatz Q , || have led to the discovery of 
a number of unusual properties (including spin-charge separation) which are fundamentally non-perturbative. 

In this paper, we analyze a quantum impurity model which can be physically realized in a resonant tunneling 
junction between multiple quantum wires or quantum Hall edges. Our interest in this problem is threefold. First, the 
results we find - both intermediate-coupling fixed points and enhanced conductance due to an analogue of Andreev 
reflection at strong-coupling - are interesting in and of themselves because they truly are, to use a cliche, exotic. 
Second, both the methods used and the result may shed light on some of the recurrent themes in the study of 
correlated electron systems in which a single-particle picture is not valid. In particular, we demonstrate a highly 
non-trivial duality that exchanges strong- and weak-coupling. Recent progress in supersymmetric field theory and 
string theory hints at the possibility that such strong-weak coupling dualties are a common, perhaps even generic, 
feature of stongly-coupled field theories. The duality discussed in this paper has a very rich structure and is one of 
the best examples of such a duality in a strongly correlated electron system. Finally, this model appears to be more 
generic and less fine-tuned than many similar ones, which leads us to hope that our findings could have consequences 
for future measurements. 

In the next section, we formulate a model describing several Luttinger liquid leads. Electrons can tunnel at a 
point-contact from one of the Luttinger liquids to a resonant state (e.g. a quantum dot or island); from the resonant 
state, they can then proceed and tunnel to another of the Luttinger liquid leads. A renormalization group analysis 
shows that when the Luttinger liquid parameter, g, is greater than 1/3, the tunneling process is relevant. In section 
III, following KSpJ], we go to a limit in which we can make an instanton gas expansion of the strong-coupling limit 
; an examination of this limit suggests a strong- weak coupling duality. This duality leads us to propose the phase 
diagram of figure 3. There are three interesting points in this phase diagram at which we can extract a more detailed 
understanding of the physics of this model. At g = 1, the electrons in the leads are non-interacting. If we assume 
that there is no interaction between the electrons at the ends of the leads and an electron on the resonant state, then 
the problem is a free fermion problem, and can be solved exactly; the solution is discussed in section IV. If, however, 
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we assume that there is such an interaction, as we do for g ^ l[], a different fixed point results. We make a conjecture 
about the relationship between these fixed points. At g — v% the model is self-dual; this property allows us to deduce 
the conductance. Finally, for g > 9, the strong-coupling fixed point is stable. At this fixed point, as we explain in 
section V, we find an analog of Andreev reflection, which leads to enhanced conductance, G > g. We also compute 
charge-transfer selection rules which elucidate the nature of this fixed point. We emphasize throughout the place of 
this model within the general framework of boundary conformal field theory and describe the most unusual features 
- namely the duality and the Andreev processes - from several different points of view. 



II. THE MODEL 



A. The Model and Formalism 



We consider a model in which N leads are coupled to each other through a resonant state, as in figure 1. One 
possible realization of this model is a quantum Hall bar in which quasiparticlcs or electrons can tunnel between 
several edges by first hopping from one edge to a dot or anti-dot and then hopping from there to another edge. An 
alternative implementation of this model is a resonant tunneling junction between TV quantum wires. The former is 
more naturally described by the 'unfolded' formalism of figure 2a in which the leads are described by chiral bosons 
on an infinite line: 

/OO /> 
dx J dT-jt- d x Md t + Vd x )<f>i (2.1) 

g and v are respectively the Luttinger paramater and velocity of the bosons, which we take, without loss of generality, 
to be the same in all leads. The quantum wire problem is more easily visualized in terms of the 'folded' setup of figue 
2b, in which the lead is modelled by a non-chiral Luttinger liquid on the half- line x < cQ 

S = f dx [ dr -!- d^d^i (2.2) 

The field, ip, is taken to have the periodicity condition^] ip = ip + 2n(2y/g). In terms of chiral fields, ip = pn + P>l- By 
this 'folding' procedure, we have essentially mapped <f>(x > 0) to ip^, as depicted in figure 2. 

The term which transfers charge to the resonant level is (in the 'unfolded' formalism, the corresponding term is the 
same, but with tpi/2^fg replaced by </>): 



■S=/ / dr^T fas+e-*^/ 2 ^ + Vj S-e^W/ 2 ^ (2.3) 



Here, we have replaced the charge state of the resonant level by a spin-I/2 degree of freedom. The spin raising and 
lowering operators, S ± , are the creation and annihilation operators for an electron or quasiparticle on the resonant 
level. The cocycles, gi, must anticommute, fa, 7?,| = 25ij, so that the tunneling operators have the correct bosonic 
commutation relations. This is true even when (|j|) transfers anyonic quasiparticles between the leads, so this model 
at g = 3 describes tunneling between quantum Hall edges via an anti-dot in the interior of a Hall droplet at v = 1/3. 
For N = 3, the rji's can be represented by Pauli matrices. When the leads arc decoupled (t — 0)', the fields ipi 
have Neumann boundary conditions at x = 0; for t ^= 0, some other conformally invariant boundary condition is 
dynamically generated in the infrared. 



1 We need such an interaction in order to pass to the Toulouse limit, as we discuss below. The g — 1 model can be continuously 
deformed into the g 7^ 1 models only when this interaction is non-vanishing. 

2 The two models are not quite equivalent since in a quantum wire or any other non-chiral Luttinger liquid, the electron 
creation operator has s pin (i.e. h — h) 1/2 and scaling dimension l/2g (i.e. h + h) while the electron creation operator in a 
chiral Luttinger liquid (2.1) is a dimension-l/2gr, spin-l/2<? operator. The two models can be mapped into each other by a 
transformation which mixes left- and right-moving modes, but point-contact tunneling is insensitive to this mixing, so all of 
our results apply equally to both the 'folded' and 'unfolded' model. 

3 The quantity g is related to the usual compactification radius r of the bosonic string via r — 2^fg. 

4 See Appendix 1. 
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Here, we are assuming that the level is perfectly resonant and that the different leads are coupled to this level with 
the same hopping strength strength, t. In an experiment, the resonance can be tuned by controlling one parameter, 
such as a backgate voltage. If there are three leads (the simplest case with a non-trivial phase diagram), then two 
more parameters must be tuned to ensure that the hopping strengths are equal. 

The fields fa, ifi can be interpreted in terms of the voltage drops along and between leads. In the 'unfolded' 
formalism, the field fa can be discontinuous across x — 0, and this discontinuity, <^(0+) — fa(0—), is proportional 
to the voltage drop across x = in the i th lead. When the leads are decoupled, there is no voltage drop along the 
'unfolded' leads, fa(0+) = fa{0—) or, equivalently, <£r(0) = <^l(0) (Neumann boundary condition). On the other 
hand, the voltage drop between leads i and j at the contact is proportional to fi(0) — <Pj(0). In most of the following, 
we will use the 'folded' formalism, but all of our results can be reinterpreted in the other language. In the appendix, 
we discuss the conventions for these bosonic fields, ip^. In particular, we discuss the mode expansions of these fields 
and the zero modes, which play a crucial role in the following analysis. In terms of the 'momentum' zero modes (see 
appendix), the Neumann boundary conditions have the effect of reflecting the zero modes of incoming states into 
those of outgoing states, P| = Pjj. When t ^ 0, these momenta are instead shifted, P[ = P l R + Q l . The allowed 
shifts, Q 1 , lie on a lattice which is connected to the problem of quantum Brownian motion in a periodic potential, as 
we will discuss in the next section. 

First, however, we note that the Kubo formula for the conductance (obtained in the usual way, see e.g. [0], by 
introducing a vector potential, A, between the resonant level and one of the leads, say lead 3, and differentiating the 
partition function with respect to A) takes the following form: 

G = 2g (l-^(<p 3 (x = 0,u)w(x = 0,-u))) (2.4) 

When t = 0, 923 is a free field (Neumann boundary condition), so|w|((^3(x = 0,uj)ip^(x = 0,u>)) = 2t: and therefore 
G = 0, as we would expect since the leads are decoupled. 



B. The Toulouse Limit and Quantum Brownian Motion 



Le t us f ocus , for the moment, on the case N = 3. The generalization to arbitrary N is straightforward. We rewrite 
and (|J) as: 



(2.5) 



where <p = (<£>i, </?2, ¥?3) and Rj 



-1,-1), R2 



1,2,-1), R3 



1, -1,2), R ± 



3V2S 



(1,1,1). 



With this notation, we have anticipated the mapping to the problem of quantum Brownian motion on a honeycomb 
lattice with lattice vectors Rjj,R|,R| and S z keeping track of the sublattice. There is one step left before such an 
identification can be complete, na mely decoupling Rj^ ■ ip by going to the Toulouse limit as in To do this, we 

modify the Hamiltonian related to (2.5) by 



H^H + t z S(x) n») S z /2^g, 



(2.6) 



where II is the momentum conjugate to (p. The added term is an interaction between the charge of the resonant level 
and the charge density of the lead at the point contact. As a result of this term, the modified Lagrangian describes 
an interacting system even at g = 1 although in this case the in tera ction takes place only at x — 0. It is natural to 



assume that this term does not affect the low-energy physics of (2.5) or, in other words, that (2.5) and the modified 
Lagrangian flow to the same infrared fixed point. This does not, appear to be the case at g = 1, but this might 
be special to g = 1. The advantage of adding such a term to the Lagrangian is that we can now |^,^| perform a 

canonical transformation ^] generated by U = e '''E i ' p '' () ' S! / 2 V9 i This has the effect of simultaneously removing the 



5 To be more rigorous, we should use U = e 1 **^* Vi(o) s^/i^/g g s *^2ji v ^ s ^/ 2 Vs ^ wmc j 1 j s overall charge neutral since only 
integer charges can be added to the systen. The second exponential compensates the fractional charge added at by removing 
an equal amount at — 00. 
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term which we just added and removing the Rj^ ■ tp terms from the exponentials in the tunneling Lagrangian if we 
choose t z — 1/N. This leaves us, finally, with the Lagrangian: 

C = i- + tS(x) (vjS + e~ m3 ]i + , h S < n k ) (2.7) 



where 



and 



[-j=(-<Pl + ipz), - + 2(^ 3 ), 4=(^1 + ¥>2 + ^3)) (2. 



Rl = yj(-^3/ 2 '-V2,0) 
R 2 = yX(V3/ 2 , -1/2,0) 



R 3 = ^i- (0,1,0) (2.9) 

Sfc^ and 9fcj, are the Cartan generators of an SC/^) that 'rotate' the leads (which is a symmetry of the free Lagrangian 
at certain special points such as g = 1, 1/2). Yi and Kane [Q showed that the 3-channel Kondo problem is one of a 
class of models (namely the g = 1/2 point) which may be formulated as the quantum Brownian motion of a particle 
on a honeycomb lattice. In ( [2.7| ), we have almost the same problem. The crucial difference is the presence of the 
rjiS which results in a 7r flux through each plaquette of the honeycomb lattice. This may be seen by considering the 
amplitude for a circuit around a plaquette, which involves the product Tjir]2r]3T}irj2T]3 = — 1. 

The RG equation for t may be obtained from the scaling dimension of the field e lR k , which is |R l | = l/3g: 

§-(•-£) « + - 

Hence, for g < 1/3, t flows to zero in the infrared and the leads are decoupled. For g > 1/3, t grows with decreasing 
energy scales. The upshot of this growth will be analyzed in the next section using a duality property of this model. 
The partition function may be expanded perturbatively in powers of t: 



Z = Y, fdn...dT n sl J2 e J R h ex P H e ^ K h ■ Ri 3 In Vi -tj\ + iir9(n - Tj ) (1 - 5 hh ) ] (2.11) 

li = 1, 2, 3, ej = ±1 and the e^'s must alternate chronologically. If we ignore the second term in the exponential, this 
is the partition function (at g = 1/2) of the 3-channel Kondo model. It is a two-component Coulomb gas. The second 
term gives a minus sign whenever the order of two unlike hops is exchanged, thereby implementing the 7r-flux. 



C. An Auxiliary Model 



We will also consider a simpler model (which is discussed in g ) for the purposes of comparison with and illumination 
of the resonant tunneling model described above. This model can be analyzed without going to a Toulouse limit, 
and it exhibits Andreev reflection at a strong coupling fixed point and a duality property with a straightforward 
interpretation. This instills us with more confidence that these properties of (2.7) are generic and are not particular 
to the Toulouse limit. It is defined by 



£ = d x <j>i{d T + id x )(j>i + tS(x) ^ 



(2.12) 



Mj = 1,2,3 



This is a model of quantum Brownian motion on a triangular lattice. In fact, this is the same triangular lattice 
which is the underlying Bravais lattice of the abov e ho neycomb lattice, a s m ay be seen by writing the Lagrangian 
( 2.12 ) as (d^k) 2 + ^(x) ^ 7 -e k ' Rj +h.c. with k as in (2.8) and Rj given by ( |2.9|) with the first and second components 
interchanged. At g = 1, (2.12) has a fermionic representation: 
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C = 4 (d T + id x ) il>i + 1 6(x) e y WjV-j (2.13) 
This is not a free fermion problem because the fc rmi on interacts with a spin-1/2 degreee of freedom, 77, which is present 



to give the correct commutation relations, as in (2.3). This model is actually a generalized multi-channel Kondo model 
in which the conduction electrons transform in an SU(2) triplet. The infrared fixed point can be solved for exactly § 
in complete analogy with the methods employed in the ordinary multi-channel Kondo model Q. Interestingly, it is 
related to the ordinary four-channel, spin 1/2 Kondo fixed point. According to ||, the model flows to strong coupling 
at g = 1. Hence, the fixed point of || is an example of the Andreev reflection phenomenon which we discuss below. 
The advantage of this model lies in the fact that there are no complications related to the Toulouse limit, as there 



are in (2.7). It is particularly simple from the point of view of duality. 



III. DUALITY 

If g > 1/3, t is a relevant coupling, so an initially small t grows in the infrared. When t is large, the interaction 

term, tS(x)J2i ^/jS ,+ e~ lRI ' k + r/jS' _ e 4R *' k ^ will be dominant and, in a semi-classical analysis, k will be localized at 

one of its minima. These minima are just the minima of the energy bands of a particle on a tight-binding honeycomb 
lattice with tt flux per plaquette. There are four such energy bands since the tt flux doubles the unit cell and since the 
honeycomb lattice, to begin with, is a triangular lattice with a two site basis. (We represent r\j by Pauli matrices Tj.) 
They correspond to the four possible S z and T3 quantum numbers. At low energies, k will be in one of the minima 
of the lowest band. These also form a honeycomb lattice; the lattice displacements - i.e. the analogs of the R;'s - on 
this honeycomb lattice are 

Ki = ^(0,l) 

K 2 = ^(V3/2,-l/2) 

K 3 = ^(-x/3/2,-l/2) (3.1) 

The partition function can be approximated by an instanton gas in which the instantons are solutions of the 
Euclidean equations of motion in which k tunnels between different minima. As usual in this class of problems [ ^| , the 
instanton gas expansion can be formulated as a Coulomb gas. There is an additional subtlety here, however: there is 
a Berry's phase associated with the instanton solutions. Details will be given in an appendix; here we merely sketch 
the derivation. Note that the minima of the lowest band surround a point at which the two lowest bands touch. The 
Berry's phase will be the same for any path surrounding this point, so we consider a path that is very close to this 

point. For such paths, ^77jS' + e _tR3 ' k + r/jS ,_ e 4RJ ' k ^ can essentially be approximated by —5k x a z — 5k y a x . Here, 

the four energy bands, acted on by rj ® S, are reduced to the two-dimensional subspace of the two lowest bands, 
acted on by the u's. 5k x , 5k y are k x , k y measured from the contact point of the two bands. As <5k traces out a path 
around 0, the spin a rotates by 27r and therefore accrues a Berry phase of 7r. Hence, the Coulomb gas defined by the 
instanton expansion is a C oulo mb gas with phases. In fact, it is of precisely the same variety as tha t de fined by the 



perturbative expansion of (2.7)! More concretely, the instanton - or strong-coupling - expansion of (2.7) is equal to 



the perturbative - or weak-coupling - expansion of 

c r> = + v S(x) £ (vjS+t K + njS < K "') (3-2) 



(Here r is the field dual to the field k of (|2.8|) in the usual way, as reviewed in Appendix C.) The v — > limit of (3.2) 



is equivalent to the t — > 00 limit of (2.7) and, conversely, the v — > 00 limit is equivalent to the t — > limit. In effect 



the duality exchanges g — > 3/g. For small v, we can obtain the RG equation for v just as we did for t above: 

*_(!-!).+... (3.3) 



dl V 9 



Combining ( |2.10 ) and ( |3.3D , we find that the t = limit is stable for g < 1/3 while the t = 00 limit is stable for 



g > 9. In the former limit, the fields k\, k^ have Neumann boundary conditions at x — (k% is decoupled, so it always 
has Neumann boundary conditions), while the r, have Dirichlet boundary conditions. In the latter limit, the situation 
is reversed. For 1/3 < g < 9, both limits are unstable and we expect a stable fixed point at intermediate coupling 
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or, in other words, a non-trivial conformally invariant boundary condition. The situation is summarized by figure 2. 
There are two intermediate coupling fixed points at which we can calculate the conductance exactly: 

(a) At g = 1, where a free fermion formulation is available for t z = 0. We do not believe that the t z = model has 
the same physics as the t z = 1/N model, but it is instructive to compare the two cases. 

(b) At g = the model is self-dual. It may be shown that the duality exchanges ^-(k x k x ) — ► 1 — ^{r x r x ) 
(and the same for k y ) as we discuss in an appendix. At the self-dual point, (k x k x ) = (r x r x ), so G = <?(2/3) = 2/V3- 
We will also discuss at length the conductance at: 

(c) The strong-coupling fixed point, which is stable for g > 9. 



First, however, we will make a few more comments on the duality between fl2.7p and (3.2). One point which should 
be emphasized is that the duality is only approximate. It is strictly a duality between the instanton gas expansion of 



(2.7) and the pert urbative expansion of (3.2) (and vice versa). In the asymptotic low-energy limit, the instanton gas 



expansion of (2.7) is the dominant contribution to the partition function when t is large, but at finite energy there 
are corrections. If we were to attempt to formulate an exact duality, thes e co rrections would be manifested by the 
presence of a presumably infinite number of additional irrelevant terms in (3.2) 

The perturbative expansion of (2.7) is an expansion in curr ent-generating charge transfer events while the instanton 
gas is an expansion in voltage-generating phase slips (in ( [3.2] ), the roles are reversed). This formulation of the duality 
concentrates on the values of the fields at the point contact. A related but alternative way of understanding this 
duality arises from the natural notion of duality inherent in th e bulk (i.e. the duality of closed strings with toroidal 
comp actifi cation) . Let us first look at the simpler model ( |2.12 ). Follow ing the same steps which led to (^^), we see 
that (2.12) is dual to a theory described by the same Lagrangian (2.12), but with the replacement g — > 3/4g [Q. Let 
us consider the finite temperature partition function of this model in a finite-size system of length L, with Neumann 
boundary condition at x — L and the interaction at x — 0, as in fig. 4a. This partition function can also be viewed (by 
turning it on its side) as the closed string amplitude for propagation between the dynamical boundary state at x = 
and the Neumann boundary state at x = L, as in fig. 4b. The closed string states are specified by their momenta, 
winding numbers, and oscillator mode occupancies (see Appendix C for a brief summary). The allowed momenta, P, 
for the fields k = (-j=(~-ipi + (p 2 ), "^( — <Pi ~ <f2 + 2953)) are determined by the condition that the operator e lP k be 
well-defined under the identification ipi = ipi + 2ir^/2g; the momenta form a triangular lattice with lattice constant 
^/l/g.[] The winding numbers, W are the set of identifications, k = k + W; they form a triangular lattice with lattice 
constant 2-^/^/3. There are two dual descriptions which result from exchanging the momenta and winding modes. 
This is precisely th e sam e duality between triangular lattices which exchanges the strong and weak coupling limits of 
(2.12). The model (2.7) can be embedded within this picture. The only additional structure is that the displacements 
on the triangular lattice (i.e. charge transfers or phase slips) are split into pairs of displacements on the honeycomb 
lattice in both the original and dual theories. Yet another interpretation in terms of S'-matrix selection rules will be 
dicussed in the context of the Dirichlet boundary condition M . 



IV. SOLUTION AT G = 1 



At g = 1, the model defined by (p.l|) and (2.3) has the free fermion representation (in particular, with t z = 0): 



C = ip\ (d T + id x ) il>i + d)d T d + t S(x) ^ + ipi$ 



(4.1) 



The creation and annihilation operators of charge on the resonant state are denoted by , d rather than , and 
{d 7 "0} = {d, ip^} = 0, {d, d^} = 1. This free fermion problem can be solved exactly. The equations of motion for ip 
and d are: 



d r ipi(x) = d x ipi(x) +tdS(x) 



(4.2) 



Integrating the first equation between x = — e and x = e and Fourier transforming, we find 



3 As usual, we ignore k z — -^(^i + <fi2 + ^3), which decouples. 



G 



ipi(w, 0+) - ipi(w, 0-) = -itd(uj) 

wd(w) = tSjil)i(uj) (4.3) 

i 

In the second equation, ipi = (^(0+) + ipi(0—))/2. From these N + 1 equations, we can extract ipi(ui,0+) and d(u;) 
in terms of ipii 1 ^, 0— ). The solution may be summarized by the 5-matrix, ipi(bj, 0+) = Sij ipj{oj, 0— ), where: 

_ JV-2 

^• = ^|, fori^j (4.4) 

The resulting conductance is 

G frccfcrmion = (JV-1) (1) (4.5) 

which, for 3 leads is G = 8/9. It is somehwat remarkable that a free fermion problem could be an interediate- 
coupling fixed point with a non-trivial conductance. However, this is the maximal possible conductance consistent 
with unitarity and permutation symmetry for a 3 lead free fermion problem. In other words, if we assume that 

= t, for i ± j (4.6) 

then unitarity, SijS*kj = 5ik, imposes the constraint |r| > 1/3, and, hence, G < 8/9. In the next section, we will 
discuss even larger conductances and the physics behind them. 

First, however, we will comment on the relationship between the t z = and t z — 1/N fixed points. We do not 
believe that they are the same for two reasons. First, we expect G/g to be non-decreasing as g is increased. While 
G(g = l,t z = 0) < G(g = V%,t z = 1/3), G(g = l,t z = 0) > G(g = Vs,t z = 1/3)/a/3. Hence, we expect that 
G(g — l,t z = 1/3) < G(g = V~3,t z = l/3)/V3 < G(g — l,t z = 0). An additional point for consideration is that a 
small t z is an irrelevant perturbation at the t z = fixed point, as may be seen by direct calculation. Similarly, a small 
deviation of t z from 1/3 is irrelevant at the t z = 1/3 fixed point, as may be shown perturbatively for g — > 1/3; it is 
reasonable to assume that this is true even at g = 1. Hence, it is plausible that the t z = fixed point described above 
lies out of the plane of the phase diagram of figure 3 with an unstable fixed point separating it from the t z = 1/3 
fixed point. 



V. DIRICHLET BOUNDARY CONDITIONS AND ANDREEV REFLECTION 

A remarkable feature of this model reveals itself when we consider the conductance at the t = oo fixed point, 



which is stable for g > 9 (and for g > 1 in the auxiliary model (2.12)). At this fixed point, k% = -^(~ ( Pi + Vi) and 
&2 = 77g( — fi ~ ¥2 + 2^3) have Dirichlet boundary conditions at x = 0, while k% = ^(fi + ¥2 + ^3) has Neumann 
boundary condition. As a result, 

H( 9 3(, = o )W Wi = o )W )) = ?M( fcl (o |W ) fcl (o )W )) + iH(yo,^ 3 (o,,)) 

= 3(0) + |(D-J ' (5.D 

where the se cond equality follows from the respective Dirichlet and Neumann boundary conditions of ki and k^. 

G™ ax = \g (5.2) 



Hence, from (^4) we have: 

4 

This is an astonishing result, since it implies that the conductance is greater than 'perfect' conductance, G — g^\ We 
interpret this as the signature of Andreev reflection: the conductance is greater than its naive maximum value because 



7 The scrupulous reader might worry that this surp risin g finding is due entirely to the Toulouse limit and is therefore incorrect. 
However, since the same conductance is found for ( |2.12| ) (which does not involve a Toulouse limit) at its strong-coupling fixed 
point, we believe that this result is robust. 
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a hole is backscattered at the point contact. Before pursuing this point further, let us note that for general N, the 
corresponding formula for the conductance at the strong-coupling fixed point is: 



GT X = 9 (2 - |) (5.3) 

For N = 2, the maximum conductance is G — g, the naive value. For N > 2, the maximum conductance is greater 
than this value, saturating at G — 2g in the TV — » 00 limit. 

Why do we say that the enhanced conductance is due to Andreev reflection? In (2.4), 2 — 2 -^(^3^3) is, essentially, 
the transmitted fraction of the incoming current; 2^-((p 3 tp 3 ) — 1 is the reflected fraction. At g — 1, where the leads 
have a free fermion description, transmission, i, and reflection, r, coefficients can be defined; 2-t^r ((^3933) — 1 = \r\ 2 . By 
charge conservation we also have (N — l)\t\ 2 + \r\ 2 = 1. G > g precisely because the reflection coefficient is negative. 
In other words, the reflected current is a negative fraction of the incoming current - i.e. it is a current of holes rather 
than electrons. 

Physically, the Dirichlct boundary condition corresponds to the limit in which there is no voltage difference between 
the different leads. For N — 1 > 1 only a fraction of the current which leaves one lead enters any one of the other 
N — 1 leads. Without Andreev scattering, this would lead to a voltage drop between the leads, but Andreev processes 
offset this voltage. An alternative perspective on the multichannel Dirichlet boundary condition is reminiscent of the 
situation explored in Q. Suppose we view N — 1 of the leads as a single, aggregate lead described by a single charge 
boson with g a ggr = g(N — 1). Then, tunneling between the remaining lead and the aggregate lead is precisely the 
problem of tunneling between dissimilar Luttinger liquids considered in || . This problem can be transformed to one 
with two identical Luttinger liquids with l/g e tt = (l/g + l/g aggr )/2 = N/2(N — l)g. For such a problem, it is not 
surprising that the maximal conductance is G max = 9eff = 2g(A — 1)/N. 

Yet another means of characterizing the Dirichlet boundary condition is by S'-matrix selection rules for soliton 
scattering at the junction. Following 0, we obtain these by rewriting the chiral fields <fii, (f>2, <^3 in terms of the dual 
fields k x , k y , which are free fields at the strong coupling Dirichlct boundary condition fixed point. Working in the 
chiral ('unfolded') notation, where k x = (<pi — (f>2)/v2, k y — (—(pi — <p2 + 2</>3)/\/6, we can define dual free fields via 

k x — k x 0(—x) — k x 6(x) 

k y = k y 9(-x) - k y 9 (x) (5.4) 



These dual fields can be identified with r x and r y which occur in the unfolded form of the dual action ( p.2|) . Now, 
we can express the in terms of the free fields k x , k v , which are just r x , r y . This allows us to calculate the matrix 
elements 



(5.5) 



The operators e Ttq j <fe( ±00 )/v / 29 create or destroy states with well-defined charges in the leads; the matrix elem ent s 



( |5.5D are proportional to the S'-matrix elements between these different charge sectors. For generic g' n ' out , (|5.5| ) 



will vanish, which means th at t here is no scattering between these charge sectors in the strong coupling (Dirichlet 



boundary condition) limit. (5J5) will be non- vanishing only if the correlation function is charge neutral for each of 



the free fields r x , r y , k z . Or, in other words, if 



q? - q 2 a = - {q° ut - q°2 ut ) 

-q? - qf + 2qf = - {-q? lt - qT" + ^T') (5-6) 



Solving for the charges of the 'out'-states, one finds that the charge transfers lie on a honeycomb lattice, 

(5.7) 



/q"'\ fq u i\ 2 

Aq= q? 



•93 




3 




Note that, for general 'in' states which carry in each lead multiples of the unit of charge, the charges of the 'out'- 
statc in the individul leads are in general no longer multiples of the unit charge. This is a phenomenon analogous to 
the N > 3 flavor Callan-Rubakov effect( p!o|). In fact, the 'auxiliary model' at g = 1, discussed at the end of Section 
II, is an example where this situation occurs at an infrared fixed point for free electron leads. For example, 
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^ e -i4>j(x=oo)/^2g e icf> 1 (-oo)/^2g \ _q 

for j — 1, 2, 3. In other words, a unit of charge cannot be transferred from one lead to another or even reflected by 
the junction! On the other hand, 

/ e i0i(a;=oo)/- v /2g e -2i02(oo)/. v /2ff e -2i03(oo)/v / 2g e 3i0i(-oo)/ v /3g \ q (5 9) 

This is a clear illustration of the Andreev reflection property of the Dirichlet boundary condition for N > 2. Three 
incoming electrons in lead I are scattered into 2 electrons into each of leads 2 and 3 and an Andreev reflected hole in 
lead 1. 
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APPENDIX A: COMMUTATION RELATIONS FOR THE KLEIN FACTORS 

We need to introduce the all-important Klein factors, rji, because we would like to treat the fields ipi as independent 
bosons which commute with each other. Since the underlying electron or quasiparticle operators are mutually fcrmionic 
or even anyonic, Klein factors must be introduced to compensate. 

We begin with the commutation relations for the chiral version of the tunneling operators, — e , which 

are obtained from those of a single chiral boson by imagining that the 3- lead dot/anti-dot setup is deformed as in 
figure 5. As a result of the chiral boson commutations relations, the tunneling operators commute since they don't 
cross (see figure 5). This holds whether the objects which tunnel arc fcrmions or anyons. 

In our model, we represent Ty by Ty = rjjS ± e^^ j where rjj and S ± commute with each other and with <fij, and 
the 4>j's are mutually commuting. To ensure that the tunneling operators commute, we must take rjirjj = —r\jT\i. 

As an aside, we note that if the tunneling paths were to cross, however, the commutation relations are modified to 
TijTki + e 2j!l / g TkiTij = 0. It is hard to imagine a setup in which this occurs, but for such a scenario, we would need 
to take the even more exotic condition 771172 = — e 2n ^ 9 ^rji an( i cyclic permutations. 



APPENDIX B: INSTANTON GAS BERRY'S PHASE CALCULATION 



As we briefly sketched in section III, when t is large, the interaction term dominates the action (2.7). If we treat 
k classically, it will be localized at one of the minima of this term. To find these minima, we need to diagonalize 

the 4x4 matrix £\ (rjjS+e-**-'* + TfjS-e^^y There are 4 solutions for each k, corresponding to the 4 bands of 

a particle on a honeycomb lattice with 7r flux per plaquette. Physically, the fourfold multiplicity is due to the two 
charge states of the resonant level and the two states of the auxiliary two-state system (i.e. 77) which keeps track of 
the statistics, while k represents the amount of charge which has been transfered between the leads. Diagonalizing, 

J2i \T)jS + e~* R% '* + r/jS^e 111 (we represent the by Pauli matrices Tj), we find the four eigenvalues: 

1/2 



e(fc) = ± \ 3 ± \l 9 - (3 + 2 cos k 2 /V2 - 2 cos(fci ^3/2.9 + k%/ y/2g) - 2 cos(fci ^3/2^ - k 2 / \/%) J (Bl) 
The minima of each of these bands form a honeycomb lattice with translation vectors 2irQi, where the Qi are given 



in (3.1) 



We now consider the instanton gas expansion of the partition function, where the instantons are solutions of the 
classical equations of motion in which k tunnels between neighboring minima. The modulus of the amplitude for 
these tunneling events can be obtained in the standard way . The phase can be obtained from the following Berry's 
phase argument. The eigenvector associated with the lowest energy band is determined by two spinors, i.e. it lies 
in the direct product space of the two two-dimensional spaces acted on by S ± and r\i. As k tunnels from minimum 
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to minimum, around a plaquette, these two spinors rotate. The phase aquired in a circuit around a plaquette is 
determined by the angles traced out by these spinors, e^ Borry = e l ^ es+e ^/ 2 , where the factor of 1/2 follows from the 
fact that S and rj are spin-1/2 degrees of freedom. Since, for any circuit, 63 and 9 V must be multiples of 27r, the only 
possible non-trivial phase is w. 

Let's consider the plaquette formed by the following six minima: (0, ±iry/2g/3), (iry/2g/3, ±27r v / 2g/3), 
(27iy2 5 /3,±7iV2<7/3). These six minima surround a maximum of the lowest band, at (%^/2g/3, 0), where the two low- 
est bands touch. The Berry phase will be the same for any loop which encloses [%^/2gJZ, 0) precisely once since such 
loops can be adiabatically deformed into each other. The Berry phase is most simply computed for an infinitesimal 
loop enclosing (7ry2g/3, 0). For such a loop, we can approximate k = (w^/2g/3, 0) + p, and 



(Vi ® S+e- 4ir - k + rjj <g> 5-e 4R! k ) « ( ~S X - ^-S y ) ® 1 + ( ^-5 X + -5„ J <8 fa. 




As p adiabatically traces out a loop enclosing p = 0, the effective Zeeman field 'seen' by ij traces out a circle but the 
effective Zeem an fi eld 'seen' by S does not. In other words 9^ — 2it while 9s — 0. This can be made more transparent 



by projecting (B2) onto the two-dimensional subspace which is degenerate at p = 0. In an orthonormal basis of the 



two eigenvectors with degenerate eigenvalues at p = 0, (|B2j) can be rewritten as: 

h{px,Py) ■ t + const. (B3) 

where r are a set of Pauli matrices and h(p x ,p v ) rotates by 2ir as p adiabatically traces out a loop enclosing p = 0. 

Hence, a tt phase is aquired in a circuit about a plaquette. Combining this with the magnit ude s of the terms in 
the standard Coulomb gas expansion for the instant on g as ffi], we see that the dual theory to (2.7) is also a theory 
defined on a honeycomb lattice with 7r-flux, namely ( |3.2|) . 



APPENDIX C: BOUNDARY CONDITIONS ON BOSONS 



In this appendix we summarize the conventions which we use for compactificd bosons. Consider a single boson 
ip(x, t) compactified on a circle of radius r. On a space of size I with periodic boundary conditions the action is 



S = i / dx I dr(d^)(d^) 



where the functional integral is to be performed under the identification 

(p(x, t) = ip(x, t) + 2irr = (p(x + I, r) 
In the Hamiltonian formalism, the field operator is 

ip(x) = <pl(x) + Vr{x) 

where 



ip L {x) ^ X L +x P L — +i2_^—=[b L . n e 1 +b' L , n e 1 \ 
1 n=i V n 

(x) = X R - x P n -j + iJ2 -j=[bR;ne- ixn2v/l + b R . n e ixn2 ^ 1 } (CI) 



VR. , 

71—1 V 



and b R]n are independent boson creation and annihilation operators, and 

[X l ,Pl] = {X R ,P R ]=i(^-) 
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are two independent zero mode coordinate and momentum operators. In these conventions, the two point funtion of 
the boson is 



< tp(x,T)ip(0,Q) >= -ln|z| 2 , (z = r + ix) 
The momentum operator, conjugate to the field (p(x) is 
2?r 



(C2) 



n(a;) = (1/2) — {(P L + Pr) + Mb L ; n e lxn27r/l + b[. n e^ n2 ^ 1 } + ]T yf^n^^ 11 + b R . n e ixn2 ^ 1 }} (C3) 



The total momentum 



n=l 



n=l 



II = / dxU(x) = (1/2) [Pi + P R \ — 
Jo » 



is canonical conjugate to the total zero mode coordinate X — Xl + Xr. Since the latter is periodic with period 2nr, 
the eigenvalues of 'dimensionless momentum' P must be of the form 



P=[P L + Pr] = 



2n 



n£ Z 



(C4) 



On the other hand, periodicity under x — > x + I gives, using (CI) 

W = Pl — Pr — rm, m G Z 



(C5) 



where we denote this quantity by the winding number W . The two conditions (C4) and (05) imply together that 



n 1 
Pl = - + -^rm; 
r 2 



n 1 
Pr = -rm: 



The (normal ordered) hamiltonian is 



H 



2n 

T 



where 



Note that 



Loa 



[Pl + Pr] 



Lq-j 



-P 2 



J2 n b R;n b ^ 

n=l 



(C6) 



(C7) 



(C8) 



[(n/r) 2 + (rm/2) 2 



gives the total scaling dimension. 

The (imaginary time evolved) Heisenberg operators are obtained from the expressions in ( Cl,p3) by 
for left movers, and by ix — > z* = r — ix for right movers. 

Of interest is also the dual boson field, 

(p{x) = (Pl(x) - <Pr(x) 
We see that the automorphism of the canonical commutation relations 

Xr ; -> -X R , Pr^-Pr, b R>n -> -b R>n b Rn ^-b Rn , r^f = 2/r 



IX —> Z = T + IX 



(C9) 



(CIO) 



(all left- movers unchanged) maps the boson field <p(x) into its dual <p(x). The dual field then satisfies the periodicity 
condition <p(x) = (p{x) + 2ixf where f is the dual compactification radius. Note that, most impo rtan tly, the duality 
automorphism exchanges the lattice of momenta P (C4) with the lattice of winding numbers W (C5). 

Boundary Conditions: Next, consider the compactified boson in semi-infinite space, < x < oo. First, we impose a 
Neumann boundary condition (b.c.) on the boson field at x = 0, 
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d x <p{x,r)^0 (x->0) (Cll) 

For the purpose of analyzing this b.c. it is convenient to view the imaginary time coordinate as a spatial coordinate 
x = r, and the original spatial coordinate x as (imaginary) time, f = — x. The new complex coordinates, are just 
rotated by 90 degrees with respect to the original ones, z = f + ix — iz, z* = f — ix — {—i)z* (this may be viewed as 
a trivial conformal transformation). We may quantize the system as before, but now on equal f slizes. The 'spatial' 
coordinate x — r is now compact, corresponding to the original system being at finite inverse temperature (3 < oo. 
The Heisenberg operators in this quantization become 

y L (x, f)=X L - i(f + *2) + i J] ^[^e^ W + fct.^-^nfe/ZJ] 



n=l 



f) = X fl -i(f-i£) P fl y +if^^[6 fl; ne (f - <£) " 2 ^ + ^.„e-^- i5 )" 2 ^] (C12) 

^ n— 1 v 

Note that the currents associated with translations of the two chiral bosons are the left/right momenta: 

2tt 2n 
Jl = idztpL = Pl^t + oscillators, Jr = idgipn = Pr.—t + oscillators 

P P 

The Neumann b.c. now becomes an identity for Heisenberg operators, acting on a boundary state \N >: 

df<p(x, f)\N (f -> 0) (C13) 

Integrating this equation over x implies in particular that the total momentum operator P = + Pr annihilates 
the Neuman n bo undary state but the winding number W = Pl — Pr = 2Pl may take on any value on the lattice 
specified in (|C5|) . (This is clear since the Neuman b.c. only constrains the derivative of the field.) Specifically, the 



only momentum states which the Neumann boundary condition supports are those with n = in (C6). Unfolding the 
semi-infinite system with a boundary thus gives an infinite chiral system of bosons with only those momenta allowed. 
Next consider the Dirichlct boundary condition 

ip(x,T)->0 (x->0) 

After 90° degree rotation, this becomes an identity for Heisenberg operators, acting on a 'Dirichlet' boundary state 
\D >: 

<p(x,f)\D>^0 (f->0) 

This implies that the total winding number operator, W — Pl — Pr annihilates the Dirichlct boundary state, but the 



total momentum P = 2Pl lies on the lattice specified in (C4). (The winding number is seen to be zero also because 



of the operator identity Xl + Xr = which follows from the Dirichlet b.c. ) S pecifically, the only momentum states 



which the Dirichlet boundary condition supports are those with m — in (C6). Unfolding the semi- infinite systems 
gives thus an infinite chiral system with only those momenta allowed. 

Duality and boundary conditions. It immediately follows from the discussion in the paragraph above that the duality 
operation exchanges Dirichlet and Neuman boundary conditions. In other words, a Dirichlet b.c. on the field ip is a 
Neuman b.c. on the dual field (p and vice versa. 

APPENDIX D: TRANSFORMATION OF THE CONDUCTANCE UNDER DUALITY 

The conductance can be obtained from the two point function 

< J 3 (t u x 1 )J 3 (t 2 ,x 2 ) > 

where 

J 3 (t,x) = JI + Jr 
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and J RL is proportional to d(p 3 R L . At the location of the point contact, we obtain four terms, which are pairwise 
equal 

< J 3 (r 1; 0)J 3 (r 2 ,0) >= 2(1 + A)/{ n - r 2 ) 2 
where A is the amplitude of the right- left current-current correlator, 

< J 3 l {t 1 ,x 1 )J 3 r {t 2 ,x 2 ) >= A/{ Zl -z* 2 f 

(zi = Ti + ixi), and 

< Jl(n,xi)Jl(T 2 ,x 2 ) >= l/4> 

(similar for right movers). 

The conductance can be written as 

G = g(l-A) 

This we may rewrite as 

G = 2. 9 [l/2 - A/2] = 2g[l - (1 + A)/2] 

On the other hand we have 

{l + A)/2= M < <p 3 {x = 0,^V 3 (x = 0,lj) > 
2tt 

where w is a real frequency, and i/? 3 = (p 3 L + ip 3 R . Under duality, A —> —A, which may also be written as 

< (p 3 {x = 0,u)ip 3 (x = 0,w) (1 - P- < cp 3 {x = 0,uj)(p 3 (x = Q,oj) >) 
2tt 2tt 

Clearly, for Neumann boundary conditions, we have have A = 1, yielding G = 0. For Dirichlet boundary condition 
on tp 3 (as opposed to Dirichlet boundary conditions on k x , k y ), A = —1, yielding, G = 2g. 

Since we will not impose boundary conditions on the fields <p J directly, we rewrite the conductance in terms of the 
fields defined in Eq(2.7), giving 

G = — [1 — — — < k x k x > aJ ] 
Under the duality transformation, k x — > k x — r Xl we have: 

O Z7T 

At the selfdual point, g = g = V3, and [1 - H < k x k x > u ] = ^ < r x r x > u = 1/2. 
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FIG. 1. A multi-lead resonant tunneling setup. 
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(a) 
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(b) 
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FIG. 2. The physical picture in the (a) 'unfolded' formalism and (b) 'folded' formalism. 
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" g ' 

FIG. 3. The phase diagram. The horizontal axis measures g, the Luttinger liquid parameter, and the vertical axis measures 
t, the hopping strength, and its dual variable, v. The RG flows are as indicated. The intermediate coupling fixed points which 
are stable for 1/3 < g < 9 are represented by the curve connecting the weak and strong-coupling fixed points at g = 1/3 and 
g — 9, respectively. 
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FIG. 4. The (a) finite-temperature partition function of this model can be represented as (b) the closed string amplitude for 
propagation between boundary states D and N. 
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